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limit. An orbital stability result for the Cauchy problem is deduced from the variational 
approach. 

Keywords: compact self-adjoint operators; trace-class operators; mixed states; occupa- 
tion numbers; Lieb-Thirring inequality; Schrodinger operator; asymptotic distribution 
of eigenvalues; free energy; temperature; entropy; Hartree-Fock model; self-consistent 
potential; orbital stability; nonlinear equation; loss of compactness 

AMS Subject Classification (2000): 35Q40 (81V45; 47G20; 81Q10; 82B10) 
1. Introduction 

The Hartree -Fo ck model with temperature is a simple extension of the Hartree- 
Fock model ^1 The minimization of the free energy determines an equilibrium 
state in the presence of a thermal noise, whose effect is to populate electronic 
states corresponding to excited energy levels. Compared to the energy functional of 
the standard Hartree-Fock model, the free energy of the Hartree-Fock model with 
temperature T has an additional term, which is the product of T with an entropy 
term. The standard Hartree-Fock model corresponds to the limi t cas e T = . 

The Hartree-Fock model with temperature has been studied I ^ I 24 | j n the case of 
an entropy based on the function (3{v) = vlogv + (1 — v) log(l — v) . The main 
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drawback in this case is that the free energy is unbounded from below when the 
model is considered in the euclidean space. 

It has been recently established that the free energy for a quantum system 
described by a mixed state or a density operator, in the presence of an external 
potential V, can be uniformly bounded from below by a functional depending on V . 
This follows from interpolation estimates of Gagliardo-Nirenberg type for systems, 
which are equivalent to Lieb-Thirring inequalities. The result of course depends on 
the convex function on which the entropy is based. We shall not make a direct use 
of such inequalities, but our approach will rely on classes of entropies for which the 
free energy is well defined and semi-bounded if V(x) = Z/\x\ . 

We shall add another ingredient, which comes from the modeling in quantum 
mechanics. We will impose that occupation numbers take values in [0, 1] because 
the electrons are fermions. This can be formalized by considering entropies based 
on convex functions with support in [0, 1] and by assuming that they take infinite 
values on R\ [0, 1] . In this paper, we shall adopt the language of density operators: a 
state of the system is represented by a trace-class self-adjoint operator 7 satisfying 
< 7 < 1 . Occupation numbers ni are the eigenvalues of the density operator: 
7 = J2i>i n i \ u i)( u i\ ■ When only a finite number of eigenvalues is non-zero, the 
operator is of finite rank. We shall say that it describes a pure state when all non- 
zero occupation numbers are equal to 1, i.e. 7 is an orthogonal projector. Minimizers 
of the Hartree-Fock model are of this form, but this will not be the case in general 
when the temperature T > is big enough: in this case the corresponding minimizer 
is a mixed state, that has some non-integer eigenvalues. 

A practical consequence of the positivity of the occupation numbers is that the 
exchange term is dominated by the Poisson term in the Hartree-Fock energy func- 
tional. In the case of a bounded domain, not much is changed when the temperature 
is turned on. A pure Hartree model (without the exchange term ], whic h is usually 
called the Schrodinger-Poisson system, has already been studied 121113] 

The purpose of this paper is to deal with the whole euclidean space and we 
shall see that many families of free energies can be considered, each of them giving 
rise to minimizers under some conditions which eventually depend on the tempera- 
ture. This is not surprising. Similar phenomena have indeed been observed in other 
models of mechanics, for instance in kinetic theory. The most studied example is 
probably the case of self-gravitating gases in astrophysics. In such models, relaxation 
mechanisms are not so well known, but stationary states can be observed and the 
corresponding free energy functional, which is also called the energy-Casimir func- 
tional, can then be recovered quite easily. Stationary solutions being characterized 
as minimizers of the free energy, orbital stability appears as a simple consequence 
for the collisionless version of the equation ^Sl. In presence of collisions, the picture 
is compatible with hydrodynamic or diffus ive limits, which provides a very natu- 
ral link with nonlinear diffusion equations ^1 Similar results are known in other 
frameworks like fluid mechanics ^ or nonlinear Schrodinger equations The strik- 
ing observation is that many different stationary states can be used simultaneously 
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for measuring the orbital stability of a given solution to the Cauchy problem. 

In this paper, we will focus on the mathematical aspects of this model. We 
consider free energy functionals of the form 

7 ^£ HF ( 7 )-TS( 7 ) 

where £ HF is the Hartree-Fock energy and the entropy takes the form 

S( 7 ) :=-tr (/%)), 

(3 being a convex function on [0, 1] . We study under which conditions (for instance 
on the total number of electrons, or global charge) such a free energy has a mini- 
mizer. 

Next we consider the stability of the associated time-dependent equation de- 
scribing the evolution of the state 7 , 

i^ = [Sy,7]> (1-1) 

which is called the von Neumann equation. Here -ff 7 is a self-adjoint operator de- 
pending on 7 but not on the function (3. Orbital stability is then a straightforward 
con sequen ce of the variational method and the existence theory for the Cauchy prob- 
lem™ Since we consider a large class of entropies based on various functions (3 , 
we thereby construct a large class of associated stable states of the time-dependent 
equation . 

The paper is organized as follows. In the next section, we state our main results 
and only give the shortest proofs. Other proofs are detailed in Section [31 



2. Results 

2.1. Definition of the free energy 

The Hartree-Fock energy, written in terms of the density matrix 7 , reads 

£ HF (7):=tr((-A) 7 )-^/ ft dx +1 D{p ^ _ II j J^l^ dx dy . 

Jrs M 2 2 JJ R 3 XR 3 \x - y\ 

(2.1) 

The first term of the right hand side of (|2.1[) is the kinetic energy of the electrons 
and the second is the electrostatic interaction with a classical nucleus of charge Z , 
located at e I 3 . We use the notation 

DU,g):=ff f^Mdxdy 

JJrzxm 3 \x-y\ 

for the classical Coulomb interaction between two densities of charge / and g . 
Such a quantity is well-defined as soon as, for instance, f,g S L 6 / 5 (IR 3 ) , by the 
Hardy-Littlewood-Sobolev inequality ^B. The last two terms of the r.h.s. of (|2.1ll 
describe the interaction between the electrons. The classical electrostatic interaction 
D(pj, pj)/2 is usually called the direct term. The last term of (|2.ip is a purely 
quantum term, which is a consequence of the Pauli exclusion principle, called the 
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exchange term. The Hartree-Fock energy (|2.ip was first studied from a m athema tical 
point of view by Lieb and Simon l2H and then by various other authors | 23 | 5 | 4 | 

The Hartree-Fock model with temperature ^ relies on the free energy, which is 
defined as 

^( 7 ):=£ HF ( 7 )+Ttr(/?( 7 )). (2.2) 

The function (3 in (|2.ip is a convex function defined on [0, 1] with values in R . To 
take into account the constraint on the occupation numbers, we extend it to +oo 
on R \ [0, 1] . The quantity — tr (/3( 7 )) is the entropy and T is the temperature. In 
quantum mechanics, a common choice is (3{v) = z^log^ + (1 — v) log(l — v) , but 
the model makes sense only on a compact subset of R 3 or for finite-rank density 
matrices ^1. We shall investigate other entropies for which the problem can be 
studied in the whole euclidean space. 

In the whole paper, we assume that 

(Al) (3 is a strictly convex C 1 function on (0, 1) , 

(A2) /3(0) = and (3 > on [0, 1] . 

Then we introduce a modified Legendre transform of (3 

g(X) := argmin (A v + (3{v)) , 

0<^<l 

i.e. 

g(X) = sup {inf{C9')- 1 (-A) ) l} , o} . 

Notice that g is a nonincreasing function with < g < 1 . Also define 

/T(A) :=Xg(X) + (f3og)(X) . 

As we will be interested in the minimization of the energy (|2.1[) under a charge 
constraint of the form tr 7 = q , we can substract T /3'(0) tr 7 to the definition of the 
free energy (|2.2p . Hence we shall assume from now on that 

(A3) (3 is a nonnegative C 1 function on [0, 1) and (3'{0) = , 

without much loss of generality. With this assumption, = and g is positive 
on (—oo, 0) . Then we also assume that (3* satisfies the following condition, 

(A4) £j 2 |/r(-Z 2 /(4Tj 2 ))| <oo. 

The role of this T-dependent condition is to ensure that the ground state free 
energy is finite, using an estimate based on the eigenvalues of the no-spin hydrogen 
atom. We recall that the eigenvalues of —A — Z/\x\ are — Z 2 /(4j 2 ) and that each 
eigenvalue is degenerated with multiplicity j 2 . 

Example 2.1. A typical example is 
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which satisfies (Al)— (A4) as long as 1 < m < 3 . In this special case 

min{(^)^ T ,l} if A<0, 

otherwise , 



.9(A) 
and 



/ 3*(A) = -(m-l)(^)" 1 - 1 if -m<A<0. 

We define next an adapted functional setting for the definition of the energy in 
terms of density matrices. Let us introduce the following space of operators 

Sj:= {7 : L 2 (M 3 ) ^ L 2 (M 3 ) | 7 = 7*, 7 G ©1 , v^A | 7 | e 61} , 



where we have denoted by 61 the space of trace-class operators 1 1 \ This is a 
Banach space when endowed with the norm 

ll7lU=tr | 7 | +tr(V^A| 7 |v^A) . 
We also introduce the convex closed subset of Sj defined by 

a: := {7 G ft I < 7 < 1} . 

Since (3 is convex and (3(0) = , we have < (3(v) < (3(1) v on [0, 1] . Hence for 
any 7 > , we have < (3(j) < (3(1) 7 . This proves that (3(j) € Si when 7661 
and therefore tr ((3(-f)) is well-defined for any 7 S K . 

For any nonnegative operator 7 , we use the shorthand notation 

tr ((-A) 7) := tr {\/~K-f eKU {+00} . 

Of course tr ((—A) 7) is finite when 7 € K . 

For any 7 € /C , we define the associated density of charge: 

p-y(x) = j(x,x) e L^M. 3 ) 

where 7(2, y) is the kernel of the trace-class operator 7 . Using the spectral decom- 
position of 7 , the following classical inequality is easily proved: 

V 7 e/C, ||VV^||™,<to(v/=A 7 V=A). 



li 2 (R3) 

Hence we have 

V 7 e/C, 11^11^1(83) < IMU • 

and, as a consequence, p 1 £ L\R 3 ) n L 3 (R 3 ) C L 6 / 5 (R 3 ). This also shows that 
Jr 3 Pj \x \ ^ ^ x ^ 00 anc ^ ^{Pi^Pi) < 00 1 by the Hardy-Littlewood-Sobolev inequal- 
ity Finally, from 7 > , it follows by the Cauchy-Schwarz inequality for se- 
quences that 

\l(x, y)\ 2 < p-y(x) p 7 (y) for a.e. (x, y) E R 3 x R 3 , (2.3) 
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SO that 

\l(x,y)\ 2 



dxdy < D(p 1 ,p 1 ) < oo . 



These considerations show that E z introduced in (|2.ip is a well-defined funtional 



on K, . We are also interested in minimizing E z under a constraint corresponding to 



a closed convex subset of K. , 

K q := {7 S K, I tr7 = q} 



2.2. The linear case 

We now recall the properties of the linear case corresponding to 

*f( 7 ):=tr((-A) 7 )-z/ £Ldx + Ttr(/J( 7 )). (2.4) 

JR3 Fl 

This corresponds to the case where the last two terms in (12.1[) are removed, ie. 

f|(7)=^(7) + ^(Ar J Ay)-| // dzdj,. 

1 1 J JWy.W F — 2/1 

n~9l 

Assume that /3 satisfies (Al)— (A4). Then a straightforward minimization 1 ^ gives 

^W = £j a /5*(^)>-oo, (2 ' 5) 

where Aj = — jjx arc the negative eigenvalues of —A — Z/\x\ with multiplicity j 2 . 
A ground state for (12. 5|) is formally given by 



7 = 5 



(t(-A-^)) 



However, this state does not necessarily belong to JC. Although its kinetic energy is 
finite by (|2.5|) . its trace 

all x (T) := tr fg (i ( - A - *,))) = £ J a ff(^) 6 (0, 00] (2.6) 

could in principle be infinite. Then the minimization problem with constraint 

mf ^(7) (2-7) 

admits a minimizcr for all q > if <7m" x = 00 , whereas it has a minimizer if and 
only if q £ [0, q l ^ x ] if gjj, n ax < 00 . In all cases, this minimizer solves the equation 



1 = 9 



(t(- a -T§i-m) 



for some fj, < , a Lagrange multiplier which is chosen to ensure that the condition 
tr 7 = q is satisfied. The function /i 1— > g(/x) := ti (g[(—A — Z/\x\ —fjt)/T)) is indeed 
non decreasing and satisfies q((i) — for /1 < 0, |/z| large enough and g(/i) — > g^ax 
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1< m < 5/3 


5/3 < m < 3 


m > 3 


&x(T) < oo 

Existence iff 
< q < gmax(T) 


gmax(T) = OO 

Existence 
Vg > 


Linear 
energy (|2.4p 
unbounded 
from 

below 


Energy (12.41) bounded from below 



Table 1. Existence and non-existence of minimizers with a finite trace for 112.71 1 for f3(i/) = u m . 

when \i — *■ . The range of the function gives all possible q's for which there is 
a minimizer for (|2.7|) . 

When gjj,ax is finite, it can take very small values, depending on the tempera- 
ture T . For instance, let us take (3{v) — v m as in Example 12.11 In this case we see 
that (j^ax < oo if and only if m < 5/3 , as summarized in TableQ] For T > Z 2 /(4 m) , 

•ten = (jfy * £ = (if^) *% ( a ^|) 

where C denotes the Riemann zeta function. We observe that q]^ ax {T) — > as 
T — > oo . We shall observe a similar result in the nonlinear case. 

2.3. Minimization of the free energy 

As a consequence of (|2.3[) and (|2.5[) , we obtain that £ z is bounded from below on K, 
for any (3 satisfying (Al)— (A4). Hence we can define 

J| (g) := inf {s§ ( 7 ) | 7 e /C and tr ( 7 ) = g} (2.9) 

and 

/f:=inf{4( 7 )| 7 e/c} = inf/' 3 ( g ). 

It is easily seen that is continuous for the topology of Sj . If we do not put any 
external potential, that is if we take Z = 0,l(j(q) can be computed explicitly: 

Lemma 2.1 (Ground state energy with Z = 0). Assume that (3 satisfies (Al)— 
(A4) for some T > 0. Then we have for any q > 

I^(q)= inf /3( ni )=/3'(0)g = 0. 

0<7li<l 

E;>i »i=9 

Proof. Let Z = . By (|2T3|) and using tr ((-A) 7 ) > , we have that 
/£(?)> inf tr (/%)) = inf = 0. 

0<n^<l — ' 
tr 7 =g Ei^xWi^a*- 1 
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Now let e > and 7 = Xa<i<_R" v i IViXv^l be a finite-rank operator such that 

tr (0(7)) < inf tr (/3( 7 ')) + e = e . 
7'eK 

tr7'=q 

As the ifi's can be chosen arbitrarily, we can assume that ifi e iJ 2 (R 3 ) for all 
1 < i < K , in which case (-A) 7 e ©1 . Let [/„ : L 2 (R 3 ) -> L 2 (JR 3 ) be the 
dilatation unitary operator defined as U v ((p)(x) := n 3 ^ 2 ip(n x) and which is such 
that U* = U 1/v . Notice that U*~fU v £ JC and that tr (f3(U* 7 U v ) = tr (/3(7)) for 
all 77 > . Using the equality (—A) C/ r) = r/ 2 [/,,(— A) , we infer 

tr((-A)E^ 7 C,) = ^ r ((- A h) ■ 

Similarly, the kernel of U*jU n is (U* -f U ri )(x,y) — v~ 3 7(x/ry, y/77) . Hence 
Pu^u^x) = 1]- 3 Pl (x/r]) and 

D(pu;-yU v ,Pu^Ur,) = ^D( Pl , Pl ) , 

// i(«;7g,)(».»)i- 1 » tMP Ws . 

77r6 7? JJrs f - 2/1 

Hence we have 

^(f/*7^)<C/^ + tr(/3( 7 )) ■ 
Taking first 77 — * 00 and then e — * yields the result. □ 

We now state our main result. 

Theorem 2.1 (Minimization for the HF model with temperature). Assume 
that (3 satisfies (Al)— (A4) for some T > . 

(1) For every q > , the following statements are equivalent: 

(i) all minimizing sequences (■y n )neN for I z (q) are precompact in K. , 
(m) I z (q) < Iz(l') f or a M 1 1 q' such that < q' < q . 

(2) Any minimizer 7 of I z (q) satisfies the self- consistent equation 

7 = 9 ((ff 7 -M)/T), fl 7 = -A-£+ Ar *|-r 1 -r^ 

R \x-y\ 

for some /1 < . 

TTie minimization problem I z {q) has no minimizer if q> 2 Z + 1 . 
(4) Problem I z always has a minimizer 7 . It satisfies the self-consistent equation 

7 = g(H,/T) . 

Remark 2.1. We shall prove below in Lemma [3731 (also see Lemma [2TTj) that 
V« , </ such that < </ < q , < I z (q - q') + I^q 1 ) = I z {q - q') . 
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The proof of Theorem 12.11 is given below in Section 13.11 Let us investigate the 
validity of Condition (ii) in Theorem 12.11 We first give a bound on the largest 
possible charge °f minimizers by comparing with the linear model. 

Proposition 2.1. Assume that (3 satisfies (Al)— (A4) for some T > 0. Then the 
minimization problem (|2.9[) has no solution if q > as defined in (|2.6[) . 

As a consequence of Proposition 12.11 and Part (3) of Theorem 12.11 we obtain 
that the largest possible charge for the nonlinear problem (|2.9j) satisfies 

fl£L<min{gi£*, 2Z+1} . 

In case of Example 12.11 for (3(v) = is m , the largest possible charge q^x converges 
to zero as T -> oo , by (|2~8|) . 



Proof. Let 7 be a minimizer of I^ z (q) for some q . By Theorem I2.1[ it solves the 
self-consistent equation 7 = g((-ff 7 — m)/^ 1 ) f° r some multiplier /i < . Next we 
notice that 

* ]— 1 - I 1 ^ ( 2 - 10 

hi \x-y\ 

in the sense of operators on L 2 (R 3 ) . In fact, we notice that 
1 ip(x) <p(y)\ 

<f(x) 2 tp(y) 2 - <p(x) <p{y) i/i(x) ip(y) 
: j ax ay > v 

\x-y\ 

by the Cauchy-Schwarz inequality. Then (|2.10p follows from the decomposition 7 = 
J2j>i n j IVj) (Vj I with Uj > . From (|2. 10|) we deduce that 

in the sense of self-adjoint operators on L 2 (R 3 ) . Since g is nonincreasing, we infer 
q = tr (<?((# 7 - M )/T)) < tr ( 5 (ff 7 /T) < tr (g((-A - Z/\x\)/T)) = ^ ax . □ 

On the other hand we can give a positive lower bound on q^ F 



imax • 



Proposition 2.2. Assume that (3 satisfies (Al)— (A4) for some T > 0. Then for 
all q such that 

0<q<min{^. 9 (^-^) , z) , (2.11) 

i>i 

Condition (ii) in Theorem \2.1\ is satisfied. 

Remark 2.2. With A^e) = -e 2 /(4j 2 ), if £\>i g{X j (e)/T) = oo for any e > , 
existence of a minimizer for Iz(q) holds for all q < Z , whereas, if Ylj>i 9(^j( £ ) /T) 
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is finite for some e > , we get the existence at least on an interval [0, (7 max ] for some 
9max > , because the right hand side of (|2.11[) does not converge to as q — > . It 
is natural to conjecture that existence holds true if we replace J2j>i di^j — o) /T) 
by J2j>i J 2 9(^j{Z ~ q)/T hi the r.h.s. of l|2.11[) . like in the linear case, but we have 
been unable to prove it. 

If (3(v) = v m with 1 < to < 3 , then Y,j>i g(Xj(e)/T) < oo for any e > , and 
so existence holds true on some interval [0, q max ] . 

Th e proof of Proposition ^. 2l is given below in Section [3~2l The case T = is well 
known I 22 | 23 | anc j ^ j s n0 ^ difg CIU t to see that estimates are uniform as T — > 0+ . 
Summarizing, we have found the following existence result. 

Corollary 2.1 (Existence of minimizers for the HF model with tempera- 
ture). Let T > . Assume that (3 satisfies (Al)— (A3), and (A4) ifT is positive. 
Then there exists q max > such that the minimization problem (|2.9|) has a mini- 
mizer for any q € [0, g max ] . 



2.4. Orbital stability 

An interesting consequence of Theorem 12.11 and Corollary 12.11 is that the set of all 
minimizers is orbitally stable for the von Neumann time-dependent equation, which 
reads 

f * 5 = i H i > > 

<^ dt (2.12) 
I 7(0) = 7o e K . 

It was proved in 013 that (|2 . 1 2[) has a global-in-time solution t j(t) S C 1 (M, K.) 
for all fixed 70 € fC q and that tr (7(4)) = q, £f F ( 7 (i)) = £f F (7o) for all * G R. 
Since (3("f) commutes with _ff 7 , it is clear that tr (/3{j{t))) is also conserved. 

For given f3, Z > 0, q > and T > 0, let M. be the set of all minimizers for 
(|2.9[) . We shall say that Ai is orbitally stable if and only if for any e > , there 
exists r) > such that for all 70 £ JC q with dist(7o,A / l) < rj , if t 1— > 7(t) is a 
solution of (|2.12p with initial data 70 , we have dist(7(t), JA) < e for all t € R . Here 
dist(7, M) := inf^e^ ||7 - 8\\sj ■ 

As a consequence of the continuity of T% and of the variational approach of 
Section l2~3l we have the 

Proposition 2.3 (Orbital stability). Assume that f3 satisfies (Al)— (A4) for 

some T > and that (i) holds true in Theorem \2.1\ for some q > . Then M. is 
orbitally stable. 

It is interesting to emphasize that we have been able to construct a large class of 
orbitaly stable states for the von Neumann equation (|2.12|) . Indeed, any function [3 
satisfying the above assumptions gives rise to an orbitaly stable set M. in /C. 
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3. Proofs of Theorem 12.11 and Proposition 12.21 

For the sake of simplicity, we shall assume that T — 1. Proving the same results for 
any T > does not add any difficulty 



3.1. Proof of Theorem [27l\ 

We recall that according to [Reed-Simon, Thm VI.26]^23, Si is the dual of the space 
of compact operators acting on L 2 (M 3 ) . Hence we can endow the Banach space Sj 
with the weak-* topology for which 7„ — 7 £ means 

tr (7„ K) -> tr (7 K) and tr (V^A-f n V^AK) -> tr (V^A 7 y/^AK) 

for all compact operators . The convex set /C is closed for the strong topology 
of io and also closed for this weak-* topology. Of course our main problem will be 
that when 7„ — 7 in K, , there could be a loss of mass at infinity in such a way that 
tr 7 < liminfn^oo tr7„ . Indeed the linear functional 7 1— > tr7 is continuous but not 
weakly-* continuous on fj , i.e. the sets 

}C q := {7 € /C I tr7 = g} 

are not closed for the weak-* topology. 

The proof of Theorem 12.11 follows some classical ideas which have been intro- 
duced in various papers H3 ^ | 3 | 17 | 18 | 

Step 1. Properties of the energy. 

Lemma 3.1. Assume that [3 satisfies (Al)—(A4) withT = 1 . Then Eg is weakly-* 
lower semi- continuous (wise-*) on IC . For every q>0, it is coercive on 

{ 7 e/C|tr7<<z}= |J JC q , . 

Q<q'<q 

Proof. Consider a sequence {7n}™eN such that j n — 7 G IC . Fatou's Lemma gives 
tr [«J— A 7 \J — A ) < liminfn^oo tr (V - A 7„ V - A ) . As {7n}neN is bounded in JC , 
7„(x, y) is bounded in i/ 1 (IR 3 x R 3 ) and yjp ln is bounded in if 1 (R 3 ) . Hence, up to 
a subsequence, we can assume that p ln — > p 7 strongly in L\ oc , 1 < p < 3 , and a.e., 
7„(x, y) — > 7(35, y) strongly in L' oc , 2 < q < 12/5 , and a.e. Hence (|2.3p and Fatou's 
Lemma give 



N - y| 



< hminf ( / / ^Wfl, (y)- | 7 n(x,y)| 2 ^ 
:R3 |x - y| 



and 

tr(/3( 7 )) <liminftr (/%„)) 
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Also 



lim / EiM dx = f BiM dx 



n— >oo 



All this shows that is wise-* on JC . Now we have 



£|(7)> Jtr((-A) 7 )+tr 

>ltr((-A) 7 )-2Z 2 



2 A bl 



which proves that f ^ is coercive as stated. □ 

Lemma 3.2. Consider a minimizing sequence {7 n }neN C K for I^il) an d assume 
that 7„ — v 7 G /C . TTien 7„ — > 7 /or the strong topology if and only if tr 7 = q . 

Proof. Assume that 7„ — 7 and tr7 = g . Then as £^ is wise-* on /C , 7 is a 
minimizer for 1% (5) . Hence linin^oo £~ (j n ) = £% (7) . By the proof of Lemma [3.U we 
see that this implies in particular linin^oo tr (V— A 7„ y/— A) = tr (\/— A 7 A ) . 
This is enough ^3 to obtain that 7n — * 7 for the strong topology of . □ 

Step 2. Binding Inequalities. 

Lemma 3.3. Assume that (3 satisfies (Al )—(A4 ) with T = 1 . For every q, q' such 
that < q' < q , we have 

4(q)<4(q'). (3.1) 

If equality holds with < q' < q , then there exists a minimizing sequence for 1% (q) 
which is not precompact. 



Proof. We consider two states 7 £ K q i and 7' £ IC q - q > such that £z(l) < IzW) + e 
-o(7') < Ifr 



and £q (7') < Iq (q — q') + e . By density of finite-rank operators in K, and of C° 



in L 2 (K 3 ) , we can assume that 7 = X^jLi n j \Pj)(<Pj\ an d 7' = XwLi n j IVjK^jl 
where the ^ and ^ are smooth functions with support in a ball 5(0, R) CM 3 . Now 
we introduce the translation unitary operator V T defined as (V T ip)(x) = ip(x — re) 
where e is a fixed vector in R 3 . We use the shorthand notation j' T := V* 7' V T . For 
t large enough, we have 77^ = 7^. 7 = , hence 7 + 7^ £ JC q and tr (/3( 7 + 7^.)) = 
tr (/?( 7 )) + tr (/9(7')) ■ An easy calculation shows that 

4(q) < 4(1 + i T ) = 4(7) + £$tf) + 0(l/r) • 

Taking first r — > 00 and then e — > yields (|3.1[) using Ig(q — q') = by Lemma |2~T1 
If there is an equality in (|3.ip for some g and < g' < q , a non compact 
minimizing sequence is built in the same way. □ 

Notice that by Assumption (A3) and Lemma |2~T| Lemma l3~3l implies that q 1— > 
Iz(q) is nonpositive, nonincreasing on [0, 00) . 
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Step 3. Proof that (i) (ii) . The implication (i) => (ii) was already proved in 
Lemma 13.31 We now prove that (ii) =4> (i) . To this end, we consider a minimizing 
sequence {j n }neti ^ K-q ■ By Lemma [3Tl {7„}„eN is bounded in ft . We can there- 
fore assume that, up to a subsequence, j n — 7 . We argue by contradiction and 
assume that j n does not strongly converge to 7 . By Lemma 13.21 this is equivalent 
to tr (7) 7^ q . Notice that by Fatou's Lemma, tr7 < lira haf n ->oo tr 7„ = q , hence 
there exists a < q' < q such that tr7 = q' . Modifying {7„} n eN if necessary, we 
can assume that {7«}neN is °f finite ran k and that (—A) 7„ € 61 for all n . 

Now we follow a truncation method ^ 7 | 18 | us choose two C°° functions x 
and £ with values in [0, 1] such that x 2 + £ 2 = 1 , x nas its support in B(0, 2) and 
X = 1 on B(0, 1) . We denote xr{%) ■= x(x/R) and £ R (x) := £{x/R) . By the IMS 
localization formula, we have 

-A = X R (-A) X R + & (-A) £ R - |V X fl| 2 - |V6?f , 
from which we deduce that 

C 

tr ((-A) 7 „) > tr ((-A) X R In Xr) + tr ((-A) £ H 7 „ - -^f ■ 
As P7- = Pxr^xr + PHr^r > we have 

D (P~t n >PlJ > £) (Pxn7 n XB>Px-R7 n Xn) + D (PiRl^R , P^rI^r ) 

where we have used that 

D (PxR7nXR=P«R7n^) = / / r — ~j d X dy > . 

j jR e \ x ~ y\ 

Using the estimates Xi?(a;) 2 ei?(y) 2 /k-y| < V# if \v\ > 3ii > and ^(v) 2 < &r(v) 2 
if \y\ < 3iZ, we obtain 

l74x,,)| 2 x^) 2 ^) 2 dxdy < i /■/■ | 7n(a , y) , 2 ^ 



l7»(^,y)| 2 Xi?.(^) 2 X3J?.(2/) 2 

We may also observe that 

l7n(a;,y)| 2 cfecft/ = tr (7 2 ) < tr (7) = q . 
The last ingredient is the 

mi 

Lemma 3.4 (Brown-Kosaki's inequality - ). Assume that (3 satisfies (Al)- 
(A3). Let 7 £ JC and consider a self- adjoint operator X : L 2 (R 3 ) — > L 2 (K 3 ) such 
that X 2 < 1 . Then 

tr(/3(X 7 X)) <tr(X0{f)X) . (3.2) 
Using (13. 2j) . we obtain 

tr (P( X R 1 XR)) + tr (/?(& 7 60) < tr (/3( 7 )) ■ 
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Hence we obtain 

£|(7n) > £z(XRlnXR.) + ^(6i7n6i) ~ HET ~ / L<\ X \< 2 R - ^ ^ " 

R<\y\<3R 1 ^' 

Notice that tr(^ fl 7„^) = q- tr (xRlnXR) and tl ' {xrIuXr}) = by 

Lemma |2. II Summarizing, we have proved that 

£>M > ^XRlnXR) - ^ - J l R<M<2R d xd y . (3.3) 

R<\y\<3R 

Lemma 3.5. We have XRln Xr ~ * XrIXR strongly in &\ , as n — ► oo . Hence 
lim tr(xi?7«Xfl) = tr(xi?7Xi?) • 

n — hdc 

Proof. We have 

tr ( X R In XR) = tr ((1 - A) 1 / 2 7n (1 - A) 1 ^ 2 (l - A)' 1 / 2 X \ (1 - A)" 1 / 2 ) , 

Since (1 - A) 1 / 2 7n (1 - A) 1 / 2 (1 - A) 1 / 2 7 (i _ A)V 2 wea kly in ft and (1 - 
A) -1 / 2 x|f (1 — A) -1 / 2 is a compact operator for any fixed R , we get 

lim tr (xRjn XR) = tr ((1 - A) 1 / 2 7 (1 - A) 1 / 2 (l - A)" 1 / 2 x 2 fl (1 - A)' 1 / 2 ) 

= tr (xfl7Xfl) ■ □ 

Passing to the limit in (13. 3p as n — * oo using the value of Jq and the weak lower 
semi-continuity of as proved in Lemma l3.ll we obtain 

R<\y\<3R 

Now we can pass to the limit as R — > oo and we obtain 

/^)>/|(g'), 

which contradicts (ii) . 

Step 4. The self-consistent equation. 

Lemma 3.6. Assume that (3 satisfies (A1)—(A4) with T = 1 and consider a 
minimizer 7 for 1% (q) . Then, for some fi < , 7 satisfies the self- consistent equation 

7 = 9{H 1 - fi) . 

Proof. The minimization problem is set on a convex set, so 7 solves the linearized 
problem 

mf {tr((tf 7 + /?'(7))7')} • 

Y£IC„ 
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The function [3 being strictly convex, 7 also solves the convex minimization problem 

inf {tr (# 7 Y + /%'))} . 
t'eK q 

This minimizer is unique and given by g{H^ — /i) where the multiplier fi is chosen 
to ensure the constraint tr 7 = q (this is an easy adaptation of (|2.5p ) . Hence 

7 = g(Hy - m) • 

Notice that er oss (i7 7 ) = [0,oo) , hence necessarily /1 < , since 7 is a trace-class 
operator and g > on (—00, 0). □ 

Step 5. Non-Existence if q > 2 Z + 1 . The proof of Lieb ^ f or the usual 
Hartree-Fock case at zero temperature applies here. We only sketch it for the con- 
venience of the reader. 

Assume that Iz(q) admits a minimizer 7 . It satisfies the self-consistent equation 
7 = g[Hry — /i) . In particular, [7, H^j = and H 1 7 < . Hence we have 

trO|-ff 7 7) < . 

Inserting the definition of H y , we obtain 

0>tr(|a;|(-A)7)-Z g + // (p 7 (x) Pl (y) - | 7 (^)| 2 ) dxdy . (3.4) 

JJe? \x-y\ 

Now we have by symmetry 



Li' 



(p-y(x) p 7 (y) - \ j{x,y)\ 2 ) dxdy 



Be f - yi 

1/7 (M + M) (p 7 ( x )p 7 ( y ) - | 7 (o;,y)| 2 ) rfsdp 

z J Jro \ x — y\ 

> \ jj ^ { Pl {x) Pl {y) - \l{x,y)\ 2 ) dxdy = l -(q 2 - tr( 7 2 )) > l -(q 2 - q) 

where we have used the fact that tr 7 2 < tr 7 = q because < 7 < 1 . Using 
that \x\ (—A) + (—A) \x\ > which is equivalent to Hardy's inequality, we obtain 
tr (|x| (-A) 7) > . Inserting in pT4jl . this yields q < 2 Z + 1 . 

Step 6. Existence for 1% . By Lemma [2~T] and Lemma l3~3l we have that q 1— > I^io) 
is non-increasing. Condition (ii) can be written Iz(q) < Iz(q') for all q' £ [0, q) . As 
we know that for q > 2 Z + 1 , I z (q) has no minimizer, necessarily is constant 
on [2Z+ l,oo): 

Vg>2Z+l, l|( ? )=/|(2Z + l). 

Hence 

l| = inf{£f(7)|7€X:, tr 7 <2Z + l} . 

Now £^ is coercive and wise on {7 G /C | tr7 < 2 Z + 1} by Lemma [3.11 hence it 
admits a minimizer, with a vanishing Lagrange multiplier fj, . 

This ends the proof of Theorem 12.11 □ 
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3.2. Proof of Proposition HOI 

We start the proof by two preliminary results. 

Lemma 3.7. For any q > , /f (<?) < and lim 9 ^ 0+ Izio) = = ^z(O) • 

Proof. By Lemma T3.31 Iz(q) is nonincreasing, so it is enough to prove the result 
for q > , small. Let us construct a state 7 £ K q with q small such that £f (7) < . 
If 95 is a normalized eigenvector of — A — Z/\x\ corresponding to the first eigenvalue 
— Z 2 /4 , we choose 

7 = q\<p){ ( p\ ■ 

The result follows from 

4(7) = + /3(g) = + o(g) as g ^ 0+ , 

since the direct and exchange term cancel for a rank-one projector. The fact that 
I§(0) = is a consequence of n g >o {7 S JC | tr (7) < q} = {0} . □ 

Lemma 3.8. Let 7 G JC q > and denote by {Xj}j>i the ordered sequence of negative 
eigenvalues of . Then we have 

{Z_q^ 2 
Tf 

In particular H 1 has infinitely many negative eigenvalues converging to zero. 



A, < . (3.5) 



Proof. The proof is inspired by several former papers I 22 | 23 [ yj e consider the sub- 
space V of L 2 (R 3 ) consisting of radial functions. For any nonnegative function 
p G V, we can estimate 

p(y) dy< ^P dx 



\x-y\ \x\ 

by Newton's Theorem. Recall that this is an easy consequence of the formula 
u(x) := / 1 dy = — pdy+ — — dy . 

J R 3 \X - y\ \x\ J\ x \<\y\ J\x\>\y\ 1 2/1 

By differentiating with respect to r = \x\, it is indeed easy to check that u is the 
unique radial solution to (r 2 u')' = —Anr 2 such that lim. r ^ ao u(r) = 0. Hence, 
denoting dX(R) the Haar measure of SO3 , if if is a radial test function, we obtain 

dxdy = / dX(R) ff vMMM dxdy 



\x-y\ Jso 3 JJrs \x-y\ 

ip 2 (x) J R3 dyj SOi dX(R) p^jRy) 
- ' \x\ 

<p 2 {x) q' , 
— n — dx 

\x\ 
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where Newton's Theorem has been applied to J go p^(Ry) d\(R) . This proves that 

iiv # 7 iiv < n y (-a - ^f) n v 

where Ily is the orthogonal projector onto V. The eigenvalues of the operator of the 
right hand side, restricted to V, are non degenerate and equal to — (Z— g') 2 /(4 j 2 ) . 
By the theory of Weyl^( we know that the essential spectrum of H 7 is [0, oo), and 
so eigenvalues can accumulate only at 0_ . □ 

Let us prove Proposition ^. 21 We argue by contradiction and consider age (0, Z] 
satisfying (|2.1ip and for which (ii) in Theorem 12.11 is false. This means that there 
exists q' € [0,q) such that I^il) = IzW) ■ By Lemma [5751 7^ is constant on [q',q] ■ 
We take q' as the smallest number satisfying this property. We notice that q' satisfies 
(|2.1ip because the r.h.s. of (|2.11[) is nondecreasing. By Lemma [5771 we have q' > . 
Also Iz(q') < Iz(q") for any q" e [0, q') . Hence q' satisfies (ii) and Iz(q') admits a 
minimizer 7 . It satisfies the self-consistent equation 

7 = 9(Hj - m) 

for some [i < . 

If n < , by Lemma [5751 we can choose an eigenvalue Xj £ ([/,, 0) with associated 
eigenfunction tp . Notice that g(Xj — p) = g'(Xj — [/,) = and 

<7 + z/|v>)(¥>| < 1 

for any v S [0, 1] . Now we compute 

f|( 7 + v = 4( 7 ) +vXj + = 4(q') +i>Xj + o{y) , 

by Assumptions (Al)— (A3). Hence Iz(q' + v) < IzW) f° r an Y v € (0, g — g') , a 
contradiction. 
If /x = 0, then 

g ' = tr( 5 (/f 7 ))>£ 5 (^f^) , 
from which we deduce that q' = q , again a contradiction. □ 
Concluding remarks 

The Hartree-Fock model with temperature is a generalization of the usual Hartree- 
Fock model, which corresponds to the zero temperature case. Many variants of this 
model can be proposed in the case of the euclidean space, each of them giving rise 
to minimizers which are in general finite rank operators, the rank being usually 
larger than the number of electrons. Generically we obtain mixed states. 

Almost for free, an orbital stability result follows from the minimization scheme. 
Proving the stability with respect to a determined stationary state is still an open 
question. This would be possible if we knew that the minimizers are isolated, a 
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property which is not granted for the zero temperature case and even known to be 
false in some special cases. 

Our framework is very natural when dealing with density operators, of finite 
or infinite rank. Recovering know results for the standard Hartrce-Fock model by 
letting the temperature go to zero is not really difficult. This should not suggest 
that the non-zero temperature case is similar to the case with zero temperature. 
For instance, the maximal number of electrons or the total charge 9^(T) that can 
be binded depends on T , and the study of this ionization threshold as a function 
of T is an open question. 
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